An Integrative Model for Soil Biogeochemistry and Methane Processes. II: Warming and Elevated CO2 Effects on Peatland CH4 Emissions

Peatlands are one of the largest natural sources for atmospheric methane (CH4), a potent greenhouse gas. Climate warming and elevated atmospheric carbon dioxide (CO2) are two important environmental factors that have been confirmed to stimulate peatland CH4 emissions; however, the mechanisms underlying enhanced emissions remain elusive. A data-model integration approach was applied to understand the CH4 processes in a northern temperate peatland under a gradient of warming and doubled atmospheric CO2 concentration. We found that warming and elevated CO2 stimulated CH4 emissions through different mechanisms. Warming initially stimulated but then suppressed vegetative productivity while stimulating soil organic matter (SOM) mineralization and dissolved organic carbon (DOC) fermentation, which led to higher acetate production and enhanced acetoclastic and hydrogenotrophic methanogenesis. Warming also enhanced surface CH4 emissions, which combined with warming-caused decreases in CH4 solubility led to slightly lower dissolved CH4 concentrations through the soil profiles. Elevated CO2 enhanced ecosystem productivity and SOM mineralization, resulting in higher DOC and acetate concentrations. Higher DOC and acetate concentrations increased acetoclastic and hydrogenotrophic methanogenesis and led to higher dissolved CH4 concentrations and CH4 emissions. Both warming and elevated CO2 had minor impacts on CH4 oxidation. A meta-analysis of warming and elevated CO2 impacts on carbon cycling in wetlands agreed well with a majority of the modeled mechanisms. This mechanistic understanding of the stimulating impacts of warming and elevated CO2 on peatland CH4 emissions enhances our predictability on the climate-ecosystem feedback

Modeling Perennial Bioenergy Crops in the E3SM Land Model (ELMv2)

Perennial bioenergy crops are increasingly important for the production of ethanol and other renewable fuels, and as part of an agricultural system that alters the climate through its impact on biogeophysical and biogeochemical properties of the terrestrial ecosystem. Few Earth System Models (ESMs) represent such crops, however. In this study, we expand the Energy Exascale Earth System Land Model to include perennial bioenergy crops with a high potential for mitigating climate change. We focus on high-productivity miscanthus and switchgrass, estimating various parameters associated with their different growth stages and performing a global sensitivity analysis to identify and optimize these parameters. The sensitivity analysis identifies five parameters associated with phenology, carbon/nitrogen allocation, stomatal conductance, and maintenance respiration as the most sensitive parameters for carbon and energy fluxes. We calibrated and validated the model against observations and found that the model closely captures the observed seasonality and the magnitude of carbon fluxes. The validated model represents the latent heat flux fairly well, but sensible heat flux for miscanthus is not well captured. Finally, we validated the model against observed leaf area index (LAI) and harvest amount and found modeled LAI captured observed seasonality, although the model underestimates LAI and harvest amount. This work provides a foundation for future ESM analyses of the interactions between perennial bioenergy crops and carbon, water, and energy dynamics in the larger Earth system, and sets the stage for studying the impact of future biofuel expansion on climate and terrestrial systems

Identification of key parameters controlling demographically structured vegetation dynamics in a land surface model: CLM4.5(FATES)

Vegetation plays an important role in regulating global carbon cycles and is a key component of the Earth system models (ESMs) that aim to project Earth\u27s future climate. In the last decade, the vegetation component within ESMs has witnessed great progress from simple “big-leaf” approaches to demographically structured approaches, which have a better representation of plant size, canopy structure, and disturbances. These demographically structured vegetation models typically have a large number of input parameters, and sensitivity analysis is needed to quantify the impact of each parameter on the model outputs for a better understanding of model behavior. In this study, we conducted a comprehensive sensitivity analysis to diagnose the Community Land Model coupled to the Functionally Assembled Terrestrial Simulator, or CLM4.5(FATES). Specifically, we quantified the first- and second-order sensitivities of the model parameters to outputs that represent simulated growth and mortality as well as carbon fluxes and stocks for a tropical site with an extent of 1×1∘. While the photosynthetic capacity parameter (Vc,max25) is found to be important for simulated carbon stocks and fluxes, we also show the importance of carbon storage and allometry parameters, which determine survival and growth strategies within the model. The parameter sensitivity changes with different sizes of trees and climate conditions. The results of this study highlight the importance of understanding the dynamics of the next generation of demographically enabled vegetation models within ESMs to improve model parameterization and structure for better model fidelity

Updated respiration routines alter spatio-temporal patterns of carbon cycling in a global land surface model

We updated the routines used to estimate leaf maintenance respiration (MR) in the Energy Land Model (ELM) using a comprehensive global respiration data base. The updated algorithm includes a temperature acclimating base rate, an updated instantaneous temperature response, and new plant functional type specific parameters. The updated MR algorithm resulted in a very large increase in global MR of 16.1 Pg (38%), but the signal was not geographically uniform. The increase was concentrated in the tropics and humid warm-temperate forests. The increase in MR led to large but proportionally smaller decreases in global net primary production (19%) and in average global leaf area index (15%). The effect on global gross primary production (GPP) was a more modest 5.7 Pg (4%). A detailed site level analysis also demonstrated a wide range of effects the updated algorithm can have on the seasonal cycle of GPP. Output from the updated and old models did not differ markedly in how closely they matched a suite of benchmarks. Given the substantial impact on the land surface carbon cycle, a neutral influence on model benchmarks, and better alignment with empirical evidence, an MR algorithm similar to the one presented here should be adopted into ELM

Evaluating Alternative Ebullition Models for Predicting Peatland Methane Emission and Its Pathways via Data–Model Fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization

Evaluating alternative ebullition models for predicting peatland methane emission and its pathways via data–model fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization

Hydrological Feedbacks on Peatland CH4 Emission Under Warming and Elevated CO2: A Modeling Study

Peatland carbon cycling is critical for the land–atmosphere exchange of greenhouse gases, particularly under changing environments. Warming and elevated atmospheric carbon dioxide (eCO2) concentrations directly enhance peatland methane (CH4) emission, and indirectly affect CH4 processes by altering hydrological conditions. An ecosystem model ELM-SPRUCE, the land model of the E3SM model, was used to understand the hydrological feedback mechanisms on CH4 emission in a temperate peatland under a warming gradient and eCO2 treatments. We found that the water table level was a critical regulator of hydrological feedbacks that affect peatland CH4 dynamics; the simulated water table levels dropped as warming intensified but slightly increased under eCO2. Evaporation and vegetation transpiration determined the water table level in peatland ecosystems. Although warming significantly stimulated CH4 emission, the hydrological feedbacks leading to a reduced water table mitigated the stimulating effects of warming on CH4 emission. The hydrological feedback for eCO2 effects was weak. The comparison between modeled results with data from a field experiment and a global synthesis of observations supports the model simulation of hydrological feedbacks in projecting CH4 flux under warming and eCO2. The ELM-SPRUCE model showed relatively small parameter-induced uncertainties on hydrological variables and their impacts on CH4 fluxes. A sensitivity analysis confirmed a strong hydrological feedback in the first three years and the feedback diminished after four years of warming. Hydrology-moderated warming impacts on CH4 cycling suggest that the indirect effect of warming on hydrological feedbacks is fundamental for accurately projecting peatland CH4 flux under climate warming

Evaluating Alternative Ebullition Models for Predicting Peatland Methane Emission and Its Pathways via Data–Model Fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial-climate feedbacks. Understanding how peatland methane (CH4) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH4 fluxes in response to changes such as elevated atmospheric CO2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH4 processes in a more mechanistic way. To address this need, we incorporated a new microbial-functional group-based CH4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH4 production—and CH4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant-mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH4 transport, respectively. A parameter sensitivity analysis showed that CH4 substrate and CH4 production were the most critical mechanisms regulating temporal patterns of surface CH4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high-carbon ecosystems from plot to regional scales

A Nonparametric Method for Separating Photosynthesis and Respiration Components in CO2 Flux Measurements

Future climate change is expected to affect ecosystem-atmosphere CO2 exchange, particularly through the influence of temperature. To date, however, few studies have shown that differences in the response of net ecosystem CO2 exchange (NEE) to temperature among ecosystems can be explained by differences in the photosynthetic and respiratory processes that compose NEE. Using a new nonparametric statistical model, we analyzed data from four forest ecosystems. We observed that differences among forests in their ability to assimilate CO2 as a function of temperature were attributable to consistent differences in the temperature dependence of photosynthesis and respiration. This observation provides empirical validation of efforts to develop models of NEE from the first-principle relationships between photosynthetic and respiratory processes and climate. Our results also showed that models of seasonal dynamics in NEE that lack specific consideration of the temperature dependence of respiration and photosynthesis are likely to carry significant uncertainties